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Abstract We consider the distribution of arg £(cr + it) on fixed lines a > A, and in 
particular the density 

d(a)= T ]ha^\{te[-T,+T] : |argC(cr + /f)| > 7t/2}\, 

and the closely related density 

d-(o)=Km^\{te[-T,+T]:3ii;(o + it)<0}\. 

Using classical results of Bohr and Jessen, we obtain an explicit expression for 
the characteristic function yfo(x) associated with arg£(tT + it). We give explicit 
expressions for d(o) and d-(o) in terms of \]f a (x)- Finally, we give a practical 
algorithm for evaluating these expressions to obtain accurate numerical values of 
d(a) and d-(a). 



1 Introduction 

Several authors, including Edwards |9] pg. 121], Gram lUTI pg. 304], Hutchin- 
son fl3l pg. 58], and Milioto J24] §2], have observed that the real part of 
the Riemann zeta-function is "usually positive". This is plausible because the 



Juan Arias de Reyna 

Universidad de Sevilla, Facultad de Matematicas, Apdo.1160, 41080-Sevilla, Spain, e-mail: 

larias@us . esl 

Richard P. Brent 

Mathematical Sciences Institute, Australian National University, Canberra, ACT 0200, Australia, 
e-mail: alf @rpb rent . com| 

Jan van de Lune 

Langebuorren 49, 9074 CH Hallum, The Netherlands, e-mail: j . vandeluneS hccnet ■ nj~| 



1 



2 



Juan Arias de Reyna, Richard P. Brent and Jan van de Lune 



Dirichlet series = 1 + 2~ 4 + 3~' v H starts with a positive term, and the other 

terms rC s may have positive or negative real part. In this paper our aim is to make 
precise the statement that 30^ (s) is "usually positive" for a := 3i(s) > \. 

Kalpokas and Steuding ifTTl . assuming the Riemann hypothesis, have given a 
sense in which the statement is also true on the critical line a — A. They showed 
that the mean value of the set of real values of £(A + it) exists and is equal to 1. 

We do not assume the Riemann hypothesis, and our results do not appear to imply 
anything about the existence or non-existence of zeros of for a > A. 

Our results depend on the classical results of Bohr and Jessen [5] concerning 
the value-distribution of in the half-plane a > A. Since Bohr and Jessen there 
have been many further results on the value distribution of various classes of L- 
functions. See, for example, Joyner |[T6l . Lamzouri |[T9l [20] |2D . Laurincikas El . 
Steuding ll27l . and Voronin OP . However, for our purposes the results of Bohr and 
Jessen are sufficient. 

After defining our notation, we summarise the relevant results of Bohr and Jessen 
in S|2] The densities d(a) and d-(o), defined in $3] can be expressed in terms of 
the characteristic function y a {x) of a certain random variable 35 associated with 
arg £(cr + it). We consider y/ a and a related function I(b,x) in Sj4]-Sf7] In TheoremQ] 
we use the results of Bohr and Jessen to obtain an explicit expression for y/o-(x). 
Theorem|2]relates logI(b,x) to certain polynomials Q„(x) which have non-negative 
integer coefficients with interesting congruence properties, and TheoremOgives an 
asymptotic expansion of I(b,x) which shows a connection between I(b,x) and the 
Bessel function Jq. Theorem|4]shows that y/ a (x) decays rapidly as x — >• °°. 

The explicit expression for y/ a is an infinite product over the primes, and con- 
verges rather slowly. In fJH]we show how the convergence can be accelerated to give 
a practical algorithm for computing y/ a (x) to high accuracy. 

In ||9]we show how d{o) and d-(o) can be computed using y/ a (x), and give the 
results of numerical computations in i fiOl Finally, in ifTTl we comment on how our 
results might be generalised. 

Elliott [ 1 ] determined the characteristic function fty (x) of a limiting distribution 
associated with a certain sequence of L-functions. We note that Elliott's *P a (x) is the 
same function as our y a (x). For a possible explanation of this coincidence, using 
the concept of analytic conductor, we refer to lfT5l Ch. 5]. Here we merely note that 
Elliott's method of proof is quite different from our proof of TheoremQ] and applies 
only to sequences of L-functions L(s,%) for which % is a non-principal Dirichlet 
character. 



Notation 

Z, Q, K, and C denote respectively the integers, rationals, reals and complex num- 
bers. The real part of z G C is denoted by 3iz, and the imaginary part by 3z. 

When considering we always have a := 3is. Unless otherwise specified, 
(7 > i is fixed. 
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Consider the open set G equal to C with cuts along (— °° + iy, j3 + iy] for each 
zero or pole j3 + iy of with j3 > A. Since is holomorphic and does not 
vanish on G, we may define log on G. We take the branch such that log is 
real and positive on ( 1 , +°° ) . On G we define arg £ (s ) by 

log CM =log|C(j)| + »-argC(i). 

P is the set of primes, and p G P is a prime. When considering a fixed prime p 
we often use the abbreviations b :~ p a and /3 := arcsin(l /b). 

\B\ or X(B) denotes the Lebesgue measure of a set B C C (orB C R). A set B C C 
is said to be Jordan-measurable if A (<9B) = 0, where dB is the boundary of B0 

2-F i (a , c; z) denotes the hypergeometric function of Gauss, see |fl] [8] . 



2 Classical results of Bohr and Jessen 

In I0J |51 Bohr and Jessen study several problems regarding the value distribution of 
the zeta function. In particular, for a > j and a given subset BcC, they consider 
the limit 

lirn ^ |{f e K : \t\ < T, log £(a + it) e B}\. 

They prove that the limit exists when B is a rectangle with sides parallel to the real 
and imaginary axes. 

Bohr and Jessen also characterize the limit. In modern terminology, they prove 
15] Erster Hauptsatz, pg. 3] the existence of a probability measure P a , absolutely 
continuous with respect to Lebesgue measure, such that for any rectangle B as above 
the limit is equal to F a (B). 

Finally, they give a description of the measure P ff . To express it in modern lan- 
guage, consider the unit circle T = {z £ C : \z\ — 1} with the usual probability mea- 
sure /i (that is d6 if we identify T with the interval [0,2n) in the usual way). Let 
P be the set of prime numbers. We may consider £2 := T p as a probability space 
with the product measure P = jj. p . Each point of £2 is a sequence co = (zp)peP, with 
each z P G T. Thus z P may be considered as a random variable. The random variables 
z p are independent and uniformly distributed on the unit circle. 

Proposition 1. Let <j > i and for each prime number q let z q be the random variable 
defined on Q such that Zq(co) = z q when CO = (zp)peP- The sum of random variables 

S:=-£log(l -/>"%)= l£^" to 4 
peP pePk=i K 

converges almost everywhere, so S is a well defined random variable. 

1 A bounded set B is Jordan-measurable if and only if for each e > we can find two finite unions 
of rectangles with sides parallel to the real and imaginary axes, say S and T, such that S C B C T 
and X(T s S) < e (see for example Halmos 1121 ). 
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Proof. The random variables Y p := — log(l — p~ a z p ) are independent. The mean 
value of each Y p is zero since 



2n 



2n Jo k 

It can be shown in a similar way that E(|F p | 2 ) ~ p~ 2 ° . Thus £ p £(|y p | 2 ) converges. 
A classical result of probability theory |[T2l Thm. B, Ch. IX] proves the convergence 
almost everywhere of the series for S. □ 

The measure P CT of Bohr and lessen is the distribution of the random variable S. 
For each Borel set B C C, we have 

F a (B) = P{(0 G £2 : S((o) G B}. 

The main result of Bohr and lessen is that, for each rectangle R with sides parallel 
to the axes, 

P ff (*) = Mm ^|{teR:|/|<r,logC(ff + ir)e*}| (D 

and the limit exists. It is easy to deduce that (Q]i is also true for each Jordan- 
measurable subset ScC, and for sets R of the form KxB, where B is a Jordan- 
measurable subset of R. 



3 Some quantities related to the argument of the zeta function 

Define a measure ji a on the Borel sets of R by fi a {B) := Pcr(R X fi). If we take a 
Jordan subset ScM, the main result of Bohr and Jessen implies that 

lia{B)= lim -L|{f eR:|/|<r,atgC(ff + i0e2?}|. 

The measure jU<j is the distribution function of the random variable 35. In fact 

jj. a (B) = P CT (R xB) = P{co G i2 : S(ffl) G R x B} = P{co G 12 : 3S(<») G B}. 
We are interested in the functions d(cr), d+(&), and d-(o) defined by 

d(a):= Hm G R : |f| < r, |argC(a + fr)| > x/2}\, 

d+(a):=jta^|{tGR:|f|<r,9lC(ff + *t)>0}| > 
rf_(a) := Hm — |{f gR: |f| < T, 9t£(CT + if) < 0}|. 
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Informally, d+{a) is the probability that 9?£(cr + /f) is positive; d-(o) is the 
probability that 3i^(a + it) is negative. We show in ClOl that d(a) is usually a 
good approximation to c/_(a). Observe that d(a) = 1 — jU(j([— jt/2,^/2]), 
d+{a) + d-(a) = 1, andd + (a) = Y.k&zVo{2kK-n/2,2kK + iz/2). 



4 The characteristic function \\r a 

Recall that the characteristic function y/(x) of a random variable Y is defined by 
the Fourier transform y(x) := E[exp(zxK)]. We omit a factor 2k in the exponent to 
agree with the statistical literature. 

Proposition 2. The characteristic function of the random variable 35 is given by 

Wc(x)=Y\l(p a ,x), (2) 

p 



I(b,x):=^- f 2 %xp(-n;arg(l-2ry' B )) dd. 



where, writing b := p° , I(b,x) is defined by 

-I 

2nJo 

Proof. By definition 

Wo(x) = J exp(«3S(o))) dco = J ]~Jexp(-LYarg(l - p~ a z P )) dm. 

By independence the integral of the product is the product of the integrals, so 

p 



a 



Each random variable z p is distributed as e' 6 on the unit circle, so 

Vato =Ili T ,I exp(-«arg(l -p-V e )) </0 = [pteV)- 



p-27T./o p 



5 The function I(b,x) 



In this section we study the function I(b,x) defined by ©. It is easy to see from © 
that I(b,x) is an even function of x. Hence, from ||2), the same is true for y a (x). 
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Proposition 3. Let b > 1 and j3 = arcsin(fe _1 ). Then 

. . If 71 ( sinf \ 

Ilb.x) = — / cos pc arctan \dt 

K Jo V b-cost J 

2b [P cos(jtf)cosf 2 f 1 / t\ dt 

-.dt = — / cos xarcsin- 



K Jo \/l-b 2 sm 2 t K Jo V bl VI -t 2 

Proof. By elementary trigonometry we find 

arg(l - b~ l e u ) = - arctan Smt — (4) 
b — cosf 

Substituting in (O gives 

, , If 271 ( sinf \ 1 [ K ( sinf \ 

Ilb.x) = — / exp /xarctan dt = — / cos \x arctan \dt. 

2k Jo V b-cost/ n Jo V b-cost/ 

To obtain the second representation, note that arctan ( sin t jib — cosf)) is increas- 
ing on the interval [0,y] and decreasing on [y,7r], where y= arccos/? -1 . We split the 
integral on [0, n] into integrals on [0, y] and [y, n\. In each of the resulting integrals 
we change variables, putting u := arctan ( sinf j(b — cosf )). Then 

f = arccos ^?sin 2 (<±cosM\/ 1 — f> 2 sin 2 M )j , 

where the sign is "+" on the first interval and "— " on the second interval. After some 
simplification, the second representation follows. The third representation follows 
by the change of variables f n- arcsin(f/f>). □ 

Lemma 1. For \t \ < 1 and all ieC, 

00 (2t) 2n "~ ' 

cos(2xarcsinf) = 2 F 1 (-x,x; \;t 2 ) = 1 + £ y -±-\\(j 2 - x 2 ). (5) 

„=1 ;=0 

Proof. In [I] eqn. 15.1.17] (also E eqn. 15.4.12]) we find the identity 

cos(2az) = iF\ {— a , a 'j 2;sin 2 z). 

Replacing a by x and z by arcsinf, we get the first half of ©. The second half follows 
from the definition of the hypergeometric function. □ 

Remark 1. An independent proof uses the fact that f(t) := cos (2x arcsinf) satisfies 
the differential equation (1 — t 2 )f"(t) — tf(t) +Ax 2 f(t) = 0, where primes denote 
differentiation with respect to f . 

Remark 2. When x £ Z, the series (O reduces to a polynomial. 
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Proposition 4. For b > 1 we have 

I(b,2x) = 2 F l (-x,x;l;b- 2 ) = 1 + £ *J\(f-f). 



-\° "■ ;=0 



Proof. From Proposition [3] we have 



7(ft,2*)=* / 1 cosf2xar C sin-^ * 



2 



x Jo V fc/ yT 

The expression of I(b,2x) as a sum follows from Lemma [T] using a well-known 
integral for the Beta function B(n + 

2 f 1 t 2n dt l , , (2n)! 



The identification of I(b,2x) as zFi(— x,x; l;b~ 2 ) then follows from the definition 
of the hypergeometric function 2F1 . □ 

Corollary 1. Ifx eZ,b 2 eQ andb > 1, then I(b,2x) G Q. 

Proof. Since I(b,2x) is even, we can assume that x > 0. Applying Euler's 
transformation 0] (15.3.4)] to the hypergeometric representation of Proposition @] 
we obtain I(b,2x) = (1 — b~ 2 ) x 2F\(— x, 1 — jc; 1; 1/(1 — b 2 )), but the series for 
2F\(— x, 1 — x; l;z) terminates, so is rational for z £ Q. □ 

We can now prove our first main result, which gives an explicit expression for 
the characteristic function y/ a defined in 



Theorem 1. For > \, the characteristic function \j/ a of Proposition\2\is the entire 
function given by the convergent infinite product 

wo{2x) = n(i + £ ± no- 2 - 2 ) • i) ■ (6) 

Proof. The identity (O follows from Proposition|2]and Proposition!?] Since Y,P 2a 
converges, the infinite product © converges for all x € C. □ 



6 The function log/(£,jt;) 

The explicit formula for y/ a given by Theorem[T]is not suitable for numerical com- 
putation because the infinite product over primes converges too slowly. In Sj8] we 
show how this difficulty can be overcome. First we need to consider the function 
logI(b,x). 
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Theorem 2. Suppose that b > max(l, |jc|). There exist even polynomials Q n (x) of 
degree 2n with Q„ (0) = and nonnegative integer coefficients q n k such that 

log/(i > 2x) = -f'^ 7 J- = -£f'^£- (7) 
By ' ' n 2 b 2 " ^fan\ 2 b 2n 

n=\ n=\k=\ 

The polynomials Q„ (x) are determined by the recurrence 



"- l 'n\/ n 



Q l (x)=x 2 , Q„ + iW = (n!)V+£^J ^. + l )Q ]+i {x)Q n -j{x). (8) 

Also, the polynomials Q n (x) satisfy 

\Qn(x)\ <«!(«- l)!max(l,|x|) 2 ". (9) 

Proof. By Proposition [4] there exist even polynomials P n with P„(0) = 0, such that 

P n (x) 1 



, n! 2 fr 2 " 



7(6,2x) = l + £ 
It follows that 

k=\ n=\ 

It is clear that expanding the powers gives a series of the desired form (|7). 

To prove the recurrence for the Qni we temporarily consider x as fixed and define 
f( y ) :=I(y- l / 2 ,2x). Then, by ®, 

iog/oo =-£%/• do) 

n=\ 

By Proposition|4]we have f(y) = 2F1 {x, —x; l;y), so f(y) satisfies the hypergeomet- 
ric differential equation 

y(i-y)f" + (i-y)f'+x 2 f = o, 

where primes denote differentiation with respect to y. Define g(y) := f'(y)/f(y)- 
Then it may be verified^ that g(y) satisfies the Riccati equation 

y(g'+g 2 )+g+- i —=o. (ID 
1 —y 

Let g(y) = YZ=aSny n , where the g n are polynomials in x, e.g. go = — x 2 . Equating 
coefficients in (fTTT i. we get the recurrence 



2 Usually a Riccati equation is reduced to a second-order linear differential equation, see for 
example Ince |14| §2.15]. We apply the standard argument in the reverse direction. 
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Y,8jgn-i-j), for »>0. (12) 

7=0 



v«+ 1 

Now, from ( fTOb and the definitions of / and g, we have 

/ \ f'(y) d d ^ Q n {x) „ 

fiy) d y dy ^ n\ l 

so we see that 



Qn 



-1 



n!(n + l)! 

Substituting $13[ in (fT2l and simplifying, we obtain the recurrence dHJ. 



(13) 



From the recurrence (|8) it is clear that Q„(x) is an even polynomial of degree 2«, 
such that 2/i (0) = 0- Writing Q„(x) = YJk=\ Qnk^K we see from the recurrence dE} 
that the coefficients q„ ^ are nonnegative integers. 

In view of ( TOT l. the inequality (0 is equivalent to |g„(x)| < max(l, |x|) 2 " +2 , 
which may be proved by induction on n, using the recurrence (TlZt . 

Finally, in view of (0, the series in (f7|i converge for b > max(l , |x|). □ 

Corollary 2. Ifb > 1, then I(b,2x) is nonzero in the disk \x\ < b. 

Proof. This follows from the convergence of the series for log I(b,2x). □ 

Proposition 5. The numbers q n £ are determined by q n \ = (n— l)\ 2 for n > 1, anJ 

n-l 

y=0 V.// Vi ' V r=JJ 



«»+!,* = L (") ( ■ " i) E 9j+l,r9n-J,k-r (14) 



/or 2 < & < n + 1, where ji = max(l,fc — n + j) and v = min(y + l,k— 1). A/so, ^„ ^ 
is a positive integer for each n > 1 anc/ 1 < k <n. 

Proof. The recurrence is obtained by equating coefficients of x 2k in (|8). Positivity 
of the <7„ j for 1 < k < n follows. □ 

Remark 3. We may consider the sum over r in (TT4l) to be over all r S Z if we define 
= for k < 1 and £ > «. The given values fi and v correspond to the nonzero 
terms of the resulting sum. 

Corollary 3. We have Yll=\ Qn,k =n\(n—\)\. 

Proof. This is easily obtained if we substitute x = 1 in the recurrence (|8). □ 
Corollary 4. We have 

q n . n = 2 2 "„!(„-l)!£4r ( 15 > 

k=\ Jo,k 

where (yo,i) is the sequence of positive zeros of the Bessel function Jq(z). 
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Table 1 The coefficients q n j,. 





1 


2 


3 


4 


5 


6 


7 


1 


1 














2 


1 


1 












3 


4 


4 


4 










4 


36 


33 


42 


33 








5 


576 


480 


648 


720 


456 






6 


14400 


10960 


14900 


18780 


17900 


9460 




7 


518400 


362880 


487200 


648240 


730800 


606480 


274800 



Proof. Define q„ : = q„,„. With k = n + l, the recurrence (TBI gives, for n>l, 

qn+i = % 0) 0+ = £ 0) - 1)**-* 1 - 

This recurrence appears in Carlitz [|7] eqn. (4)], where it is shown that the solution 
satisfies lfl5V □ 

Remark 4. The sequence (q„) is A002190 in Sloane's on-line encyclopedia of 
integer sequences (OEIS), where the generating function — log (Jo(2s/x)) is given. 
The numbers q„ enjoy remarkable congruence properties. In fact, (1131 is analogous 
to Euler's identity |Z?2„| = 2(2n)!£^° =1 (27i:£)~ 2 ", and the numbers q„ are analogous 
to Bernoulli numbers. We refer to Carlitz |]7] for further discussion. 

Remark 5. There are other recurrences giving the polynomials Q„ and the numbers 
q n k. We omit discussion of them here due to space limitations. 



7 Bounds and asymptotic expansions 

Since I(b,x) is an even function of x, there is no loss of generality in assuming 
that x > when giving bounds or asymptotic results for I(b,x). This simplifies the 
statement of the results. Similarly remarks apply to which is also an even 

function. 

Consider the first representation of I(b,x) in Proposition^ If b is large, then 

/ sin0 \ sin0 , 

arctan = — + {b' 2 ) . 

\b — cos9 J b v ' 

However, it is well-known ll32l §2.2] that the Bessel function Jq(x) has an integral 
representation 

1 f n 

Jo(x) = — / cos(xsm9)d0 . (16) 
it Jo 
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Thus, we expect I(b,x) to be approximated in some sense by Jo(x/b). A more de- 
tailed analysis confirms this (see Proposition Q and Corollary |9). The connection 
with Bessel functions makes Corollary |4]less surprising than it first appears. 

Proposition 6. For allb>\ and x£M.,we have \I(b,x) \ < 1. 

Proof. This follows from the final integral in Proposition |3] □ 

Lemma 2. For t e [0, 1] and c\ = %jl—\< 0.5708, we have 

< arcsin(/) — t < c\t 3 . 

Proof. Let f(t) = (arcsin(r) — t) jt 3 . We see from the Taylor series that f(t) is non- 
negative and increasing in [0, 1]. Thus sup fG [ n /(f) = /(l) = k/2 — 1. □ 

Lemma 3. Suppose b > 1, t € [0, 1], and c\ as in Lemma\2} Then 

< foarcsin(*/fo) -* < c^/b 2 . 

Proof. Replace f by t/b in Lemma [2] and multiply both sides of the resulting 
inequality by b. □ 

Proposition 7. Suppose b> I, x>0, and c 2 = (2 -A/%) /3 < 0.2423. Then 

\I{b 7 x)-J Q (x/b)\<c 2 x/b 3 . 

Proof. From the last integral of Proposition [3] we have 

2 f l ( t\ dt 

I{b,bx) = — / cos I oxarcsin - 



n Jq V />/ Vl -f 2 

Also, from the integral representation (TTrJi for 7o, we see that 

2 , , dt 

Jo(x) = — cos(xf) 
?r Jo 

Thus, by subtraction, 



I(b,bx)-Jo(x) = l f l f(b,x,t) 
% Jo 



dt 



VT 



(17) 



where f(b,x,t) — cos(bxarcsin(t /b)) — cos(xf ). Using |cosx — cosy)| < \x — y\, we 
have 

\f(b,x,t)\ < \bxaicsm(t/b)—xt\. 

Thus, from Lemma[3] 

\f(b,x,t)\<c lt 3 x/b 2 . 

Taking norms in (fTTI i gives 
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i, ( ,, fa) _ W i<l£/;_^£_. <18 > 

The integral in (TT~8T > is easily seen to have the value 2/3. Thus, replacing x by x/b 
in (fT~ST > completes the proof. □ 

Corollary 5. Ifb >l,x>0,C2 as in Proposition^ and C3 = \/2/n < 0.7979, then 
\I(b,x)\<c 2 x/b 3 +c 3 (b/x)^ 2 . (19) 

Proo/ It is known (TJ 9.2.28-9.2.31] that |7 (x)| < ^2j{nx) for real, positive x. 
Thus, the result follows from Proposition [7] □ 

Remark 6. The crossover point in Corollary |5]is for b w x 3 /7 : the first term in (fT9t 
dominates if <C x 3 / 7 ; the second term dominates if b x 3 ' 7 . 

Corollary 6. 7/x > 1 and b > x 1 / 2 , then 

\I(b,x)\<c 3 (b/x) l / 2 (l+c 5 b- 1/2 ), 

where c% C3 are as above, and C5 = C2/C3 < 0.3037. 

Proof. From Corollary [5] we have 

|/(MI < ca(*A) 1/2 (i +c 5 x 3 / 2 /^ 7/2 )- 

The condition b > x l l 2 implies that x 3 / 2 /b 7 / 2 < b~ l l 2 . □ 

For the remainder of this section we write j3 := arcsin(l /b). 
Proposition 8. For b > 1 a«c/ rea/ positive x, we have 

, , w r „ VF^Tcosh M -/sinh M . 

7(6,jc) = -91 ( — — / e~ A " tfa I . (20) 

' 1 — (coshM + Wb 2 — 1 sinhw) 2 



% Jo 



Proof. From the second integral in Proposition [3] we get I(b,x) = 3U(b,x), where 

2b [P irt cosf 
J(b,x) = — e M - =dt. 
71 Jo Vl-^sin 2 ? 

The function 1 — fr 2 sin 2 ? has zeros at t = ±j3 +kn with teZ and only at these 
points. Also, [5 = arcsin(l/£>) € (Q,n/2). Hence, if £2 denotes the complex plane 
C with two cuts along the half-lines (— j3] and [/J,+°°), then the function 

(cos t ) I v 1 — b 2 sin 2 f is analytic on . We consider the branch that is real and pos- 
itive in the interval (0,/3). We apply Cauchy's Theorem to the half strip 3f > 0, 
< % < j3, obtaining 
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i f°° _ xl< 2£cosh M ie ix P f°° _ m 2bcos(fi+iu) 
J(b,x) = - / e x " - du / e x " -.du. 

"Jo Vl+b 2 sinh 2 u * Jo Jl - & sin 2 (fi+iu) 



The first integral does not contribute to the real part. Taking the real part of the 
second integral and simplifying gives J20t . □ 

In the following theorem we give an asymptotic expansion of I(b,x). 

Theorem 3. For b > 1 fixed and real x — > +°° there is an asymptotic expansion of 
I(b,x). If fi = arcsin(l /b), the first three terms are given by 

2 (b 2 -\) l l A 

/(M= 

(b 2 + 2)(b 2 -l)- 1 / 4 . . _ /A . 
4V2tt x i ' 2 

(% 4 -28fc 2 +4) (b 2 -l)- 3 / 4 to lA . f 1 

64\/27r x 5 / 2 y H 11 \xV 2 



Proof. We apply the Laplace method and Watson's Lemma 1251 Ch. 3, pg. 71] to 
the representation (|20t . □ 

Corollary 7. For fixed b > 1, the function I{b,x) has infinitely many real zeros. 

Proof. This is immediate from the first term of the asymptotic expansion above. The 

M 

4 



zeros are near the points ± + kit) / ft for k £ Z>q. □ 



Corollary 8. For fixed <7 > A, the function yfa(x) has infinitely many real zeros. 
Proof. This is immediate from Proposition|2]and Corollary 
Corollary 9. 7f£> > 1 and j3 = arcsin(l /£>), then for real x +°° we have 

i(b, x ) = p l ' 2 {b 2 - i) 1/4 /o03*) + e»(^- 3 / 2 ). 

Proof. The Bessel function /oM nas an asymptotic expansion which gives 

/o(x)= (—\ (cos{x~it/4) + ^-sm{x-it/A) + o(\ 
\7tx J \ 8x \x z 

Therefore, from TheoremQ the difference I(b,x) — fi l l 2 (b 2 — l)'/ 4 ./o(j3x) is of the 
order indicated. □ 

Now we give a bound on the function I(b,x) which is sharper than Corollary |5]in 
the region x ^> b 1 ^. 

Proposition 9. For b > y/2 and real x > 5, we have 

\I(b,x)\ < I.l5l2y/b/x. (21) 
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Proof. We consider the representation (120V Take A := y/b 2 — 1 so the condition 
> s/2 implies that A > 1. It can be shown that, for A > 1 and real u > 0, the 
inequality 

A cosh m — / sinh u 



yl — (coshM + ;A sinhw) 2 



c 4 \/A 



l( M l/2,l) 



holds. Here, the optimal constant is C4 = ^/coth(2) < 1.0185, attained atA = i/= 1. 
(We omit details of the proof, which is elementary but tedious.) Hence, from d20l> . 

\I(b,x)\ < 7 ^-{b 2 -\yl A {j\- l l 2 e- xu du + J™ e- xu du} 
2C4VB 



< — - — \y%/x + e x /x 

Forx>5 we have (2c 4 /7t)(y /r 7t/x + e' x /x) < 1.1512/0:. □ 

Remark 7. The constant 1.1512 in d2"TT i can be reduced if we do not ask for unifor- 
mity in b. From Theorem[3] we have 



as x — > +°°, 



so the constant can be reduced to C3 = (2/k) x I 2 < 0.7979 for all x > xo(b). 

The following conjecture is consistent with our analytic results, for example 
Corollary |5] and Theorem [3] and with extensive numerical evidence. 

[2b 

Conjecture 1. For all b > 1 andx > 0, we have |/(fo,x)| < \ 

V Kx 

To conclude this section, we give a bound on y/ a (x). 

Theorem 4. Let o > \ be fixed. Then \y a {x)\ < \for allxG K. Also, there exists a 
positive constant c > 0.47 and xq((j) such that 

( cx l /° \ 
|V<tWI <exp -- — —j— forallreal x>x Q {a). 
y log(x l ' a ) J 

Proof. The first inequality is immediate from the definition of y/ CT (x) as the charac- 
teristic function of a random variable. 

To prove the last inequality, it is convenient to write y := x x l a . Let S^iy) be the 
set of primes p in the interval (y'/ 2 ,y]. We can assume that yf a {x) 7^ 0, because 
otherwise the inequality is trivial. From Proposition |6]and Corollary [6] we have 

Wo{x)\< n i'(p ff .*)i< n (c3( P /y) a/2 (i+c5P- a/2 )), 

which implies 
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-log|vr«,(*)|> D (-log(c 3 ) + f (logy- log/,)) +0(y l -°l 2 ). 

Using log(c3) < —0.22 and a > 1/2 gives 
-log| V /aW|>(7r(y)-7r(y I / 2 ))(flogy + 0.22)-f £ log/, + C?(y 3 / 4 ), (22) 

where, as usual, 7t(y) denotes the number of primes in the interval [1 ,y] . 
From standard results on the distribution of primes [28], we have 

n{y) = -^ + ^^ + o(^ T -] and V logp =y + O ( ) . 
logy log 2 y UgV pe % y) \\o£y) 

Substituting in (1221) . we see that the leading terms of order y cancel, leaving 

logy Vlog 2 y/ 

Since S + 0.22 > 0.47, the Theorem follows, provided y is sufficiently large. □ 

Remark 8. We find numerically that, for a £ (0.5, 1.1), we can take c = 1 andxo = 5 
in Theorem[4] 



8 An algorithm for computing \\f G (x) 

There is a well-known technique, going back at least to Wrench |33l , for accurately 
computing certain sums/products over primes. The idea is to express what we want 
to compute in terms of the prime zeta function 

*W :=!>"' (*(')>!)■ 
p 

The prime zeta function can be computed from log using Mobius inversion: 

P^) = L— log CM- (23) 

r=l r 

In fact, d23b gives the analytic continuation of P(s) in the half-plane 3is > (see 
Titchmarsh |29j §9.5]), but we only need to compute P(s) for real s > 1. 

To illustrate the technique, temporarily ignore questions of convergence. From 
Theorem |2] we have 
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io g /(p v) =-i ^i P -^. 

n= 1 

Thus, taking logarithms in ©, 

log wo(x) = -I i 9nW± p -2n* = £ Q^0_ P{2na) , (24) 

p n=\ n - n=l n - 

Unfortunately, this approach fails, because y/ a (x) has (infinitely many) real zeros - 
see Corollary |8] In fact, the series (l24l converges for |x| < |xi(ct)|, where xi(cr) is 
the zero of y a (x) closest to the origin, and diverges for \x\ > \x\ (cr)|. 

Fortunately, a simple modification of the approach avoids this difficulty. Instead 
of considering a product over all primes, we consider the product over sufficiently 
large primes, say p > pa(x, a). Corollary |2]guarantees that I(p a ,x) has no zeros in 
the disk |x| < 2p° . Thus, to evaluate yf a (x) for given a and x, we should choose 
2p° > |x|, that is po > |jc/2| 1 / cr . In practice, to ensure rapid convergence, we might 
choose po somewhat larger, say po « |4x| 1 / cr . 

For the primes p < po, we avoid logarithms and compute I(p a ,x) directly from 
the hypergeometric series of Proposition 2] 

To summarize, the algorithm for computing y/ a (x) with absolute error O(e), for 
x £ K, is as follows. 



Algorithm for the characteristic function (x) 

1. P0 ^n^i 1/a i- 

/ N { n-\ \ 

2. A i- ] [ 1 + 2na \2 il^ 2 _ ( x / 2 ) 2 ) » where is sufficiently large that 

P<po V "=l P H - j=0 J 
the error in truncating the sum is O(e). [Here A is the product over primes < pq.] 

3. B <- exp ( - Qn(yX { 2 ^ { P{2na) - £ p' 2na \ J , where N' is sufficiently 

V «=1 { P<P0 )J 

large that the error in truncating the sum is 0(e), and Q n (x/2) is evaluated using 
the recurrence (|8). [Here B is the product over primes > po.] 

4. return Ax B. 
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Remarks on the algorithm for y a (x) 

1. At step [5J P(2no) can be evaluated using equation d23k time can be saved by 
precomputing the required values £(rs). 

2. It is assumed that the computation is performed in floating-point arithmetic with 
sufficiently high precision and exponent range J6] Ch. 3], For efficiency the pre- 
cision should be varied dynamically as required, for example, to compensate for 
cancellation when summing the hypergeometric series at step 12 or when com- 
puting the term {P(2no) - Y,p< Po p~ 2 " a } at step[3j 

3. At step[3jan alternative is to evaluate Q„(x/2) using a table of coefficients q n 
these can be computed in advance using the recurrence of Proposition This 
saves time (especially if many evaluations of \j/ a (x) at different points x are re- 
quired, as is the case when evaluating d(a)), at the expense of space and the 
requirement to estimate N' in advance. 

4. The algorithm runs in polynomial time, in the sense that the number of bit- 
operations required to compute y/o(x) with absolute error 0(e) is bounded by 
a polynomial (depending on a andx) in log(l/e). 



9 Evaluation of J (a) and d-(c) 

In this section we show how the densities d(a) and d-(a) of fJJJcan be expressed 
in terms of the characteristic function y a . 

Proposition 10. For o > 1, the support of the measure )i a of$3\is contained in the 
compact interval [—L(a),L(a)] t where 



Proof. Recall that jj, a is the distribution of the random variable 35 considered in i|2] 
From (0), 35 is equal to the sum of terms — arctan((sinf)/(/? CT — cosf )) whose 
values are contained in the interval [— arcsin/?~ a ,arcsin/;>~ CT ]. Therefore the range 



Remark 9. It maybe shown that the support of fi a is exactly the interval [L(a),L(a)]. 

Remark 10. As in van de Lune l23l . we define do t° be the (unique) real root in 
(l,+oo) of the equation L(a) = It/2, and Ci to be the real root in (1,+°°) of 
L(g) = 37T/2. These constants are relevant in ClOl 



arcsin(p a ). 



p 




□ 
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Proposition 11. For a > \, 

d (a) = \-l /Vtosinf?)-. (25) 



K JO 

Proof. Recall from ^3]that 1 — d(a) = ^i a ([-n/2,n/2\). Since \jf a is the charac- 
teristic function associated to the distribution [i a , a standard resulj^] in probability 
theory gives 

1 f x exp(ixn/2) - exp(-a7r/2) 
1 — <i(cr) = — - hm / — — — — — w a (x)dx. 

Since yf a (x) is an even function, we obtain (f25t . □ 



To evaluate d(a) numerically from (125) . we have to perform a numerical integra- 
tion. The following theorem shows that the integral may be replaced by a rapidly- 
converging sum if a > I. 



Theorem 5. Let a > 1 and £ > max(n/2,L(a)). Then we have 

r 2 



K 2 ^ 1 (Kn\ . (Ml 1 
d(a) = l ) -Wa\ — sin 

2£ n ^ n \ £ J \2£ 



(26) 



Proof. Consider the function p (x) equal to p a (x) in the interval [—£, £] . Now extend 
p(x) to the real line K, making it periodic with period 2£. Thus 

P(x) = L /» ex P {^~J~) ' where f" = ^ exp ( _ ~T~) dx ' 

Now p(x) = Po(x) for |;t| < £ and p a (x) = for |x| > £. Therefore 

f - = hL po{x)exv (-ir) <fa= 27/ R p - (x)exp (-^r) ^tMt) ■ 

Since y/<j (x) is an even function, 

P W = S I >( T ) «P (—J = ^ + J L H-) co- — ■ (27) 



Now d(ff) = l-^ CT ([-7r/2,7r/2]) = 1 - f*£, 2 p a {t)dt. Since 7r/2 < £, we may 
replace p a (t) by p(f) in the integral. Hence, multiplying the equality (|27| | by the 
characteristic function of [— n/2, ti/2] and integrating, we get (l26l l. □ 



3 Attributed to Paul Levy. 
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Remark 11. The sum in d26l l can be seen as a numerical quadrature to approximate 
the integral in (|25l l, taking a Riemann sum with stepsize h = However, we 
emphasise that d26l i is exact under the conditions stated in Theorem This is a 
consequence of the measure [i c having finite support when a > l.Ifce (j , 1] then 
jj, a no longer has finite support and (|26| | only gives an approximation; however, 
this approximation converges rapidly to the exact result as I — s- °o, because ji a is 
well-approximated by measures with finite support. 

Remark 12. If we take m :=Aij% in the Theorem|5] we get the slightly simpler form 



d{a) = 1 - 




for m > max(2,M(a)), where M{a) = 4L(a)/n. A good choice if L(a) < n is 
m = 4; then only the odd terms in the sum d28l l contribute. 



Computation of c/ (a) 

Recall that c/-(cj) is the probability that 3i£((J + it) < 0. Let a^ = ak(o) be the 
probability that |arg£(<7 + /f)| > (2k+ \)%/2, that is 

a k :=\-ii a ({-(k+\)K,{k+\)n\). 

Then 

d-{a) = £(fl2t-fl2W-i) = £(-l)*fl*- (29) 

We have seen that, for a > 1 and m > max(2,4L(<r) / n), eqn. d28l gives «o = d(a). 
Similarly, under the same conditions we have 

, Ak + 2 2 " 1 (An\ . ({4k + 2)nn\ 

a k = \ V-Yo — sin -i ' — . (30) 

m 7T « \ni / \ m J 

Using < f29b and ( f30b in conjuction with an algorithm for the computation of we 
can compute d-(a) and also, of course, d+{o) = 1 —d-(o). If <7 G (i,l] then we 
can take the limit of d30b as m — > °°, or use an analogue of PropositionQ~T] to evaluate 
the constants a^. 
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Table 2 d[a) for various a 6 (0.5, 1.165] 



a 


d(a) 


0.5+ 10-" 


0.6533592249148917497 




0.5+10- 5 


0.4962734204446697434 




0.6 


7.9202919267432753125 x 10 


-2 


0.7 


2.5228782796068962969 x 10 


-2 


0.8 


5.1401888600187247641 x 10 


-3 


0.9 


3.1401743610642112427 x 10 


-4 


1.0 


3.7886623606688718671 x 10 


-7 


1.1 


6.3088749952505014038 x 10" 


-22 


1.15 


1.3815328080907034247 x 10" 


103 


1.16 


1.1172074815779368125 x 10" 


194 


1.165 


1.2798207752318534603 x 10" 


283 



10 Numerical results 

In O we described a computation of the first fifty intervals (t > 0) on which 
+ it) takes negative values. The first such interval occurs for t w 682112.9, 
and has length ps 0.05. From the lengths of the first fifty intervals we estimated that 
d-{\) pa 3.85 x 10~ 7 . We also mentioned a Monte Carlo computation which gave 
d-{\) w 3.80 x 10~ 7 . The correct value is 3.7886... x 10~ 7 . The difficulty of im- 
proving the accuracy of these computations or of extending them to other values of 
a was one motivation for the analytic approach of the present paper. 

The algorithm of Sj8]was implemented independently by two of us, using in one 
case Mathematica and in the other Magma. The Mathematica implementation pre- 
computes a table of coefficients q„ the Magma implementation uses the recurrence 
for the polynomials Q„ directly. The results obtained by both implementations are in 
agreement, and also agree (up to the expected statistical error) with results obtained 
by the Monte Carlo method in the region 0.6 < a < 1.1 where the latter method is 
feasible. 

Table 2 gives some computed values of d(a) for a £ (0.5, 1.165]. From van de 
Lune 123 we know that d{a) = d-(o) = for a > a ~ 1.19234. Table 2 shows 
that d{a) is very small for a close to ob- F° r example, d(a) < 10~ 100 for a > 1.15. 
The small size of d(o) makes the computation difficult for a > 1.15. We need to 
compute i//(j (An / m) to more than 1 00 decimal places to compensate for cancellation 
in the sum d28b . in order to get any significant figures in d(a). 

Selberg (26) (see also EOS ED) showed that, for t - unif(T,2T), 

fagC(l/2 + *) , v { „ (31) 

^logiogr 

as T —> oo, withX,T ^N(Q, 1). This implies that d(l/2) = 1, but gives no indication 
of the speed of convergence of d(o) as a J, i. Table 2 shows that convergence is 
very slow - for a — J > 10 -11 we have d(o) < |. 
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Table 3 The difference d(a) - d-{a) 



a 


d{p)-d-{p) 


0.5 + KT 11 
0.6 
0.7 
0.8 


0.1547533823 
8.073328981 x 10-" 
2.676004882 x lO" 32 
7.655052120 x lO" 210 



It appears from numerical computations that 9t£(l/2 + it) is "usually positive" 
for those values of t for which computation is feasible. This is illustrated by several 
of the Figures in J2J. Because the function ^/loglogT grows so slowly, the region 
that is feasible for computation may not show the typical behaviour of + it) for 
large t on or close to the critical line a = \ . 

Table 3 gives the difference d(a) — d-(o). For a > 0.8, there is no appre- 
ciable difference between d(a) and c/_((7). This is because the probability that 
|arg£((7 + /f)| > 37r/2 is very small in this region. Indeed, d(a) = d-(a) for all 
a > <J\ ps 1.0068, where <3\ is the positive real root of L(a) = 3n/2. 

There is an appreciable difference between d{a) and d-{a) very close to the 
critical line. For example, rf_(0.5 + 10" 11 ) ps 0.4986058426, but d{Q.5 + lO" 11 ) ps 
0.6533592249. Our numerical results suggest that lim^i nd- (o) = 1 /2. 

It is plausible that d-(j) = d+(j) = j, but Selberg's result (I3H does not seem 
to be strong enough to imply this. 



11 Conclusion 

We have shown a precise sense in which 'KC(s) is "usually positive" in the half- 
plane a — > i, given an explicit expression for the characteristic function y/a, 
and given a feasible algorithm for the accurate computation of y a , and consequently 
for the computation of the densities d(p) and d-(o). 

Our results could be generalised to cover Dirichlet L-functions because the char- 
acter X(p) m me Euler product 

p 

can be absorbed into the random variable z P whenever = 1. Thus, it would 

only be necessary to omit, from sums/products over primes, all primes p for which 
X(p) is zero, i.e. the finite number of primes that divide the modulus of the L- 
function. This would, of course, change the numerical results. Nevertheless, we 
expect 3iL(s,x) to be "usually positive" for 9?(s) > j. 
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